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The dependence of tlie tliree lowest order spatial correlation functions of a harmonically con- 
fined Bose gas on temperature and interaction strength is presented at equilibrium. Our analysis is 
based on a stochastic Langevin equation for the order parameter of a weakly-interacting gas. Com- 
parison of the predicted first order correlation functions to those of appropriate mean field theories 
demonstrates the potentially crucial role of density fiuctuations on the equilibrium coherence length. 
Furthermore, the change in both coherence length and shape of the correlation function, from gaus- 
sian to exponential, with increasing temperature is quantified. Moreover, the presented results for 
higher order correlation functions are shown to be in agreeement with existing predictions. Ap- 
propriate consideration of density-density correlations is shown to facilitate a precise determination 
of quasi-condensate density profiles, providing an alternative approach to the bimodal density fits 
typically used experimentally. 
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I. INTRODUCTION 

The appearance of coherence in a system, which is 
central to our understanding of laser and matter wave 
physics, depends on the interplay between various pa- 
rameters, such as temperature, interaction strength, con- 
finement, and dimensionahty. In particular, in the re- 
cently achieved confined weakly-interacting quasi-one- 
dimensional Bose gases in harmonic traps 
dipole tra p s and atom chips |^ llollll , 
[la . [l7l llM coherence can only be maintained across 
the entire spatial extent of the system at a sufficiently 
low temperature, which is however much lower than the 
'critical temperature' for the onset of quantum degener- 
acy. At intermediate temperatures, long wavelength fluc- 
tuations in the ph ase, restrict the coherence to smaller 

regions [13 |M S 111 111 IM 111 111 111 IM IM 113 

This effect has already been observed experimentally in 
very elongated three-dimensional Bose-Einstein Conden- 
sates ^A,^3.,M,M!,2M- 

In his seminal work , Glauber characterised the co- 
herence of a system by means of a set of normalised cor- 
relation functions. Typically, the few 'lowest order' cor- 
relation functions are enough to characterise the coher- 
ence of a system [sl Isl] . Perhaps the most important 
correlation function is the first order correlation func- 
tion, referred to in condensed matter literature as the 
off-diagonal one-body density matrix. In an atomic gas 
above the degeneracy temperature, this quantity decays 
rapidly to zero, on a scale comparable to the atomic de 
Broglie wavelength (see, e.g., |3l|). In the opposite ex- 
treme, when this quantity tends to a nonzero plateau, 
the system becomes coherent and is said to contain 
a Bose-Einstein Condensate (BEC). In degenerate one- 
dimensional (ID) gases, however, an intermediate state 
exists, in which the off-diagonal one-body density matrix 
does decay to zero within the system size, but at a much 
slower rate. In this regime, the system is said to contain 
a quasi-condensate |23 . 



The phase coherence properties of such trapped quasi- 
condensates were first investigated by Petrov et al. |21| in 
the low temperature limit, where quasi-condensate deple- 
tion can be treated as negligible. A number of alternative 
mean field approaches have appeared in the literature 
since 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 40, 4lJ. 
In particular, the theory of Andersen et al. (2J] treats 
both phase and density fluctuations in a self-consistent 
manner, thus providing the natural extension of the the- 
ory of Petrov et al. [21| to finite temperatures. Using 
this theory, we constructed the universal one-dimensional 
(ID) phase diagram in the weakly-interacting regime [23| . 
We also showed that inclusion of density fluctuations can 
potentially lead to a significant decrease in the equilib- 
rium coherence length of a trapped quasi-condensate [25j | . 
Although initial experiments performed to address this 
question found no such effect within their experimental 
resolution H, HjQ , we argued that the apparent in- 
sensitivity of the (suitably scaled) coherence length on 
density fluctuations only arises in the regime <C Tc, 
in which all experiments to date have been performed, 
where is the characteristic temperature for phase fluc- 
tuations, and Tc is the 'critical temperature' for the onset 
of (quasi)condensation. In fact, despite a relatively weak 
coupling between density and phase fluctuations in most 
experiments, the first evidence of such coupling may have 
actually been recently observed |34|. The optimum con- 
ditions for the unequivocal experimental demonstration 
of such coupling between density and phase fluctuations 
was laid out in [23 . 

In this paper, we present results for the three low- 
est order spatial correlation functions at equilibrium, 
based on the stochastic Langevin treatment of Stoof 
mii^. Where appropriate, these results are compared 
and contrasted to predictions of various mean field the- 
ories mmillliHlil. Firstly, we discuss the depen- 
dence of a suitably-defined equilibrium coherence length 
on temperature. Predictions of the present treatment 
are shown to be in full agreement with the mean field 
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theory of Andersen et al. |22|, which includes density 
fluctuations in an ab initio manner, thus complement- 
ing and extending earlier related work Further- 
more, the present work provides additional evidence to 
support our earlier claim (2^ (based on mean field the- 
ory) that the coherence length is indeed sensitive to den- 
sity fluctuations and quasi-condensate depletion in re- 
alistic parameter regimes. The spatial variation of the 
off-diagonal first order correlation function with distance 
z from the trap centre (z = 0) is evaluated in two differ- 
ent ways, corresponding experimentally to the particular 
atomic interference experiment that could be performed: 
In the first approach, atoms at the centre are numeri- 
cally 'interfered' with atoms at a distance z from the trap 
centre, whereas the other approach 'interferes' atoms lo- 
cated symmetrically about the trap centre, i.e., at points 
±z/2. Good agreement is found between these two ap- 
proaches within the quasi-condensate region. The inter- 
play between quasi-condensation and 'true' condensation 
is known to be manifested, not only in the change of the 
scaled coherence length, but also in the shape of the first 
order correlation function, which changes smoothly from 
exponential to gaussian [29|. Analysing the correlation 
functions obtained by the stochastic approach at different 
temperatures, we model the interplay between such be- 
haviour by two parameters, the first one controlling the 
shape of the function, and the second representing an 
appropriately-defined coherence length, which decreases 
uniformly with increasing temperature. 

Higher order density-density correlation functions con- 
tain crucial information regarding the system's coher- 
ence, and determination of such correlations (indirectly 
via measurement of the collisional interaction energy 
E3,E1)E3> or the three-body inelastic loss rate [soUsUl 
has played a key role in unequivocally demonstrating the 
experimental observation of BEG in ultracold 3D bosonic 
atomic gases. Moreover, second order correlation func- 
tions have been well-studied in the context of photon 
statistics in the Hanbury-Brown-Twiss experiment |52l |. 
to discriminate between thermal and coherent light fields 
[37L Is^ l. and a similar technique was recently used to 
experimentally determine the phase coherence length of 
a trapped elongated quasi-condensate 0, . Further- 
more, spatial two-body correlations have been observed 
in expanding atom clouds |53l Is^ and atoms produced 
from molecular breakup |55'|, whereas temporal correla- 
tions were measured in continuously outcouplcd atom 
lasers j5j|. This paper additionally investigates in de- 
tail the dependence of the second and third order spa- 
tial correlation functions in confined weakly-interacting 
ID Bose gases on temperature and interaction strength. 
The second order correlation function has already been 
discussed in the weakly-interacting regime in [43, with 
the crossover to the strongly-interacting regime presented 
in Is^l , and our results reveal good agreement with 
such treatments in the appropriate regimes. 

Density-density correlations are also used to directly 
determine the amount of quasi-condensation present in 



a system at any given instant. Although already ap- 
plied to study quasi-condensate growth on an atom chip 
[53 , the details of this technique were not explicitly pre- 
sented in our earlier work. In our opinion, this approach 
provides the most direct method that can be experimen- 
tally applied to obtain an unequivocal determination of 
quasi-condensate density profiles. In addition to provid- 
ing an alternative to the conventional bimodal density 
fits, the advantage of experimentally obtaining the quasi- 
condensate density by the proposed method is that it is 
not subject to any potential limitations of the particu- 
lar theory used to analyze the experimental data. This 
technique can be readily applied to existing experiments 
which performed in situ measurements of density fluctu- 
ations of quasi- ID Bose gases at equilibrium |35l| . 

This paper is structured as follows: Sec. II briefly out- 
lines the stochastic approach and other mean field theo- 
ries to which the results are compared (see also Appendix 
A) . Sec. Ill discusses the first order correlation function, 
focusing on the dependence of an appropriately defined 
coherence length on temperature (Sec. Ill A), and on 
method of evaluation (Sec. Ill B); moreover. Sec. Ill C 
analyses the crossover in the shape of the first order cor- 
relation function with decreasing temperature. Sec. IV 
discusses the dependence of higher-order density-density 
correlation functions on position and temperature (Sec. 

IV A) and effective interaction strength (Sec. IV B). Sec. 

V expounds how measurement of density-density correla- 
tions can provide an unequivocal determination of quasi- 
condensate density, which is demonstrated by means of 
suitable examples. Finally, Sec. VI summarises the main 
results of this paper. 



II. LOW-DIMENSIONAL THEORIES 

We start by reviewing the stochastic approach which 
will be used to evaluate all correlation functions pre- 
sented in this work, and also briefiy outlining the par- 
ticular mean field theories against which our results will 
be tested. 



A. Stochastic Langevin Approach 



The stochastic Langevin theory of Stoof |43l l44j is a 
non-equilibrium approach, in which the system dynamics 
is obtained via a Langevin equation governing the evo- 
lution of the quasi-condensate order parameter <I>(z,i), 
given by 0| 



9$(z,t) 



2to 



+ y'='^'(z) -^-ii?(z,i) 



^iz,t) + r,{z,t) 



(1) 



In the present work, this equation is applied to study the 
growth of a ID degenerate atomic gas in contact with 
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a thermal cloud which acts as its heat bath. The nu- 
merical implementation of this scheme, along with fur- 
ther details, are discussed in Refs. 0, while al- 
ternative but related approaches can be found in Refs. 
61 62,MJ- 

In addition to the usual kinetic, potential and interac- 
tion terms appearing within the square brackets of Eq. 
(1), Stoof's Langevin equation contains the contribution 
iR{z, t), which describes pumping of the one-dimensional 
gas from the surrounding thermal reservoir. In the clas- 
sical approximation imposed here for numerical simplic- 
ity, whereby the non-condensed ID atomic cloud relaxes 
to the 'classical' value N{£) = [f3{e — , the above 
pumping term obeys 



iR{z,t) = 



2m 



+ V'^''\z)-^i + g\<^{z,t)\■ 



(2) 



The term 77(2, t) corresponds to associated noise arising 
from the random nature of collisions occuring in the sys- 
tem. This noise term has Gaussian correlations of the 
form 



{r^*{z,t)ri{z\t')) 



i:^{z)5{z- z')5{t-t') , (3) 



where (...) denotes averaging over the realizations of the 
noise. The quantities of Eqs. (2) and (3) depend on 
the one-dimensional Keldysh self-energy, HYj^ {z), which 
accounts both for collisions that send a thermal atom 
into the quasi-condensate, and those which promote a 
quasi-condensate atom into the thermal cloud. The for- 
mulation of this theory ensures that the trapped gas 
relaxes to the correct equilibrium, in accordance with 
the fluctuation-dissipation theorem. The quantity fx 
appearing in Eqs. (l)-(2) describes the chemical po- 
tential of the system, whereas the external potential 
yext^^^ is treated as harmonic throughout this work, i.e. 
yext^^^ _ ma;^z^/2, where the corresponding confin- 
ing frequency. 

The quantity $(2,^) contains information about both 
mean field and fluctuations around it. It therefore im- 
plicitly includes both density and phase fluctuations, 
and acts as an approximation to the Bose field opera- 
tor ^{z,t). Predictions of this theory are based upon 
averaging the Langevin field $(z, t) over different noise 
realizations. In particular, this theory enables an ac- 
curate determination of both diagonal and off-diagonal, 
spatial and temporal correlation functions of any order 
at any spatial and temporal coordinate. This is achieved 
by suitable numerical autocorrelation measurements of 
the 'order parameter' <I>(z, i), upon making the following 
identification 



{^\z,t)i>{z' .t')) ($*(z,t)$(z',t')) 



(4) 



In this paper we restrict our analysis to spatial corre- 
lations at equilibrium. To study these, we first evolve the 



system for a sufBciently long time, t^q, such that it relaxes 
to the correct equilibrium, before performing the desired 
same-time autocorrelation measurements. The first order 
normalised off-diagonal correlation function g^^^{0, z) as- 
sociated with fluctuations in the quasi-condensate phase 
is thus calculated with respect to the centre of the trap 



Higher order correlation functions are also routinely ob- 
tained by a simple generalisation of the above formula. 
In particular, the second and third order correlation func- 
tions evaluated at the same spatial coordinate at equilib- 
rium are numerically obtained via 



(6) 



and 



5^^) (z) = 9^'^ {z,z,z,z,z,z, Uq) = ^^[I'llfP^^ ^ (7) 
respectively. 

The above technique is very powerful, as it additionally 
enables a full non-equilibrium determination of coherence 
properties of the system, as demonstrated in |59j |. This 
is of direct relevance to recent growth experiments 
which will be investigated in future work. 



B. Mean Field Theories 

In our initial analysis, results of the first order cor- 
relation function, evaluated from Eq. (5), will be com- 
pared against correponding results of a modified finite 
temperature mean field theory valid in low dimensions 
[2^ . Although good agreement between these two the- 
ories has already been demonstrated |23|, our aim here 
is to show that the agreement is sufficiently sensitive to 
quantitatively assess the role of density fluctuations on 
equilibrium coherence properties. In particular, under 
suitable realistic conditions, predictions of the stochastic 
approach are found to be in full agreement with the mean 
field theory which includes density fiuctuations . Im- 
portantly, however, the same stochastic predictions are 
shown to differ substantially from the limiting case of 
that mean field theory, in which density fluctuations are 
ignored This is an important remark, given that 

all coherence experiments (except l63]) have so far been 
analyzed in terms of theories which ignore density fluc- 
tuations. 

Both mean field theories mentioned above are sum- 
marised in Appendix A. For our present purposes, it is 
sufficient to quote here the final expression for the first- 
order normalized correlation function g^^\0, z), which is 



4 



given in terms of correlations of the phase operator, x{z) 
by 

gW(0,z) = e-^<K(^)-*(")l'> . (8) 

The exponent of the above expression can be written in 
ID as an appropriate sum over Legendre polynomials, 
with an additional prefactor which explicitly depends 
on the spatial extent of the quasi-condensate (see Eq. 
(Al)). The latter parameter is appropriately defined by a 
'temperature-dependent Thomas- Fermi radius', i?TF(2^), 
whose size decreases with increasing temperature. 

The low temperature theory of Petrov et al. 0| is 
obtained as a limiting case of this theory upon setting 
the quasicondensate depiction to zero, as demonstrated 
in [25| . To perform a calculation which is consistent with 
the 'classical' Langevin theory employed here, we again 
allow numerically the thermal part to relax to the classi- 
cal value N{£) = [j3{e — 

Let us now discuss the results of the theories presented 
above. 



III. FIRST ORDER CORRELATION FUNCTION 
AT TRAP CENTRE 

A. Temperature Dependence of g''-^\0,z) 

The spatial correlation function .g'^-'(0, z) evaluated at 
the trap centre at equilibrium by the stochastic approach 
is shown for two different temperatures by the brown 
(grey) lines in Fig. l(a)-(b). Corresponding predictions 
based on the theory of Andersen et al. [22|, i.e. in- 
cluding densiy fluctuations, and the theory of Petrov et 
al. |2lj |. i.e. excluding density fluctuations, are shown 
respectively by the dashed and solid black lines. The 
spatial extent of all plotted curves is restricted to points 
within the temperature-dependent Thomas-Fermi (TF) 
radius i?TF(T'), which characterises the system size at 
each temperature (see also 13)- expected, the 
correlation function in the absence of density fluctua- 
tions is consistently higher than both other predictions 
(stochastic theory and mean field theory with density 
fluctuations), with this discrepancy becoming significant 
at higher temperatures. Importantly, these latter two 
theories show consistent behaviour over the entire tem- 
perature range, with only slight differences in the shapes 
of their predicted correlation functions at intermediate 
temperatures, presumably arising because the stochastic 
treatment does not predict gaussian behaviour close to 
the origin. 

The good agreement between the stochastic approach 
and the mean field theory with density fluctuations be- 
comes more evident in Figs. l(c)-(d), which depict the 
temperature dependence of the correlation function suit- 
ably characterised by means of two independent determi- 
nations. Throughout this work, temperature is scaled to 
the 'phase coherence' temperature, = N{fuvzy / n [2ll |. 



which marks the onset of phase fluctuations and roughly 
separates the regions of quasi-condensation and 'true' 
condensation. Here N denotes the quasi-condensate 
atom number, uj^ the conflning frequency, and /i the ID 
chemical potential. 

The flrst approach we use here to compare the predic- 
tions for this correlation function over the entire temper- 
ature range, based on these three theories, is to investi- 
gate how the value of the correlation function halfway to 
the edge of the quasi-condensate, i.e. at z = Rtf{T)/2, 
varies as a function of temperature. This is shown in 
Fig. 1(c) against scaled temperature T/T^. In the 
quasi-condensate regime considered here, this quantity 
decreases smoothly to zero at a rate which depends crit- 
ically on whether density fluctuations are included (cir- 
cles, open squares), or not (filled squares), with both 
theories including density fluctuations showing very good 
agreement. 

An alternative method for characterizing the temper- 
ature dependence of the correlation function is to appro- 
priately extract a coherence length from it, scale this to 
the system size, and study its variation as a function of 
scaled temperature T/T^. A graph of this form is of- 
ten compiled by experimentalists |34j . Here we choose 
to deflne the coherence length, Lcoh, as the value of 
z at which the correlation function decays to half its 
original value, i.e. g^^HOj^coh) — 0.5. We shall deal 
with a more suitable deflnition of the coherence length 
which also accounts for the changing shape of the cor- 
relation function in Sec. Ill C. In our numerical re- 
sults, Lcoh is scaled to the zero-temperature system size 
i?TF(0), and its dependence on scaled temperature T/T^ 
is shown in Fig. 1(d). Importantly, we flnd that all 
curves which include density fluctuations lie on a univer- 
sal curve Lcoh/ Rtf{0) = exp(— T/T^), whereas curves 
without density fluctuations (flUed squares) lie consider- 
ably higher, indicating an overestimate of the actual equi- 
librium coherence length. As argued elsewhere '255, the 
shift between curves which include the coupling of phase 
and density fluctuations to those that do not, depends 
critically on the ratio of the phase coherence temperature 
T^, to the ID 'transition temperature' Tc- Importantly, 
however, we note that the trend of the observed decrease 
of the coherence length due to density fluctuations is con- 
sistent with recent experimental flndings [3^ . 

All results presented in this paper are plotted in terms 
of suitably scaled parameters (lengths, temperatures), 
such that any dependence on the particular system de- 
tails is removed. The underlying simulations were per- 
formed at constant chemical potential /i = SOftw^ for 
a ID gas of ^^Na atoms, of scattering length — 
2.75nm, under longitudinal conflnement of frequency 
ujz — "^TT X 3.5 Hz. The ID interaction strength, de- 
noted by g in Eq. (1), is given hy g — Ai:h^K,/m, 
where n = (rn/2TTh)LU j^a^D is the ID 'coupling constant' 
(or 'inverse scattering length') [b^. In our ID simula- 
tions the value of k is fixed by our choice of transverse 
confinement, taken here to be harmonic with frequency 



5 





1 




0.8 




0.6 






0.4 










0.2 







(c) 


1 




0.8 






p 


0.6 


tr" 


0.4 


o 


0.2 














4 ^ T = 0.5 T, 

Rtf(T) ~ Rtf(O) 




FIG. 1: (color online) (a)-(b) Normalized first order spatial 
correlation function g^^\0,z) at equilibrium, obtained from 
the stochastic approach (solid brown/grey lines), and from 
mean-field theory with (dashed black) and without (solid 
black) inclusion of density fluctuations. Plotted temperatures 
correspond to T/T^ « (a) 0.5, and (b) 2.8, with position z 
scaled to the effective temperature-dependent Thomas-Fermi 
radius i?TF(r), the value of which decreases from (a) to (b). 
(The corresponding density profiles are shown in Fig. 6(a)- 
(b).) (c) Dependence of (;'^'(0, z = Rtf(T) /2) on scaled tem- 
perature T/T^ based on the stochastic approach (open circles) 
and modified mean field theory with (open squares) and with- 
out (filled squares) density fluctuations, (d) Dependence of 
the scaled coherence length Lcoh/ Rtf{0) on scaled tempera- 
ture T/Tcj, for the 3 theories shown in (c). Here Lcoh is defined 
by S'''^'(0, Lco]i/Rtf{0)) ~ 0.5, where /?tf(0) is the spatial ex- 
tent of the condensate at T = 0. Only points corresponding 
to temperatures T > are plotted, since at lower tem- 

peratures the correlation function remains consistently above 
the value of 0.5, indicating full coherence. Green (grey) dia- 
monds correspond to the curve Lcoh/i?Tp(0) = exp{—T/T^). 
Plotted lines in (c)-(d) connect data points as a guide to the 
eye. 



IV ± = 2tt X 120Hz. Performing the simulation at fixed 
chemical potential leads to a slight variation in the atom 
numbers in the trap, in the range 18, 500 — 24, 000, for 
the considered temperature range, lOnK < T < 400 nK. 
An analogous variation would also be expected in experi- 
mental realisations at different temperatures and approx- 
imately fixed atom numbers. Accordingly, in the con- 
sidered temperature range, the phase coherence temper- 
ature increases approximately linearly with increas- 
ing atom number (or, equivalently here, temperature T), 
from its initial value ^ 100 nK to approximately 
Ta ~ 125 nK. 



B. Symmetrically Evaluated Correlation Function 

g^-'\-z/2,z/2) 

A number of current experiments are better suited 
to determining the correlation function symmetrically 
about a point, as done in the Hannover |3^ and Orsay 
experiments Our next task is therefore to investi- 
gate the extent to which the predictions discussed above 
depend on the precise method by which the correlation 
function is obtained. We will do this by comparing cor- 
responding predictions at the trap center, based on the 
stochastic theory. 

The 'symmetric' correlation function at the trap centre 
is thus obtained at equilibrium via the following numer- 
ical autocorrelation 



{<i>*{~z/2,t,,Mz/2,t,g)) 
V(|$(-^/2,te,)P>(|$(^/2,<e,)|^ 



(9) 

which should be compared and contrasted to Eq. (5). 

Correlation functions computed symmetrically via Eq. 
(9) are shown by the brown (grey) lines in Fig. 2(a) for 
four different values of T/T^p. These are compared to 
corresponding correlation functions ^^^•'(O, z) of Eq. (5) 
(black). We find that, for a given temperature which fixes 
the density profile, these two correlation functions arc in- 
distinguishable in a broad region close to the trap centre. 
Some differences do, however, arise towards the edge of 
the quasi-condensate, and beyond, for certain tempera- 
tures just below the phase coherence temperature T^. A 
slower decrease of the symmetrically computed correla- 
tion function (j^^-* (— z/2, z/2) is anticipated in this region, 
since the different methods of evaluation used imply that 
the considered correlation functions will decay to zero at 
different points, respectively at z « i?TF(T) for g'^^\0, z) 
and z « 2i?TF(T) for 5(1) (-z/2, z/2). 

For sufficiently low temperatures, as in Fig, 2(a)(i), 
the agreement between the two different determinations 
of the correlation function is accurate almost up to the 
edges of the quasi-condensate (i.e., z = i?TF(r)), and 
the same holds deeply in the quasi-condensate regime 
shown in Fig. 2(a)(iv), for which the two independently 
determined correlation functions are practically indistin- 
guishable. In fact, this equivalence remains nearly perfect 
over the entire temperature range, for spatial coordinates 
z < i?TF(r)/2. This is evident in Fig. 2(b) which com- 
pares the temperature dependence of the value of the 
correlation function at z = Rtf{T)/2, based on the two 
different determinations mentioned above. The reader 
is reminded that a similar determination of the coher- 
ence length was used in Fig. 1(c) to demonstrate the im- 
portance of density fluctuations. For completeness, Fig. 
2(c) plots the effective system size, Rtf{T), versus scaled 
temperature T/T^, demonstrating an approximately lin- 
ear decrease with increasing temperature. 

One might also consider performing the alternative 
analysis of Fig. 1(d), whereby the coherence length is de- 
fined by the spatial coordinate z, at which the first order 
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correlation function 5'^' decays to the value of 0.5, with 
this value scaled to i?TF(0)- In this case, at low tempera- 
tures T <T^, the two correlation functions are sensitive 
to their precise evaluation method, and the agreement 
between them would not be as good. 
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FIG. 2: (color online) (a) (i)-(iv) Spatial dependence of the 
normalised first order correlation function, calculated via the 
stochastic theory, under different autocorrelation measure- 
ments. Black lines denote the function (j'^-'(0,z) of Eq. (5), 
whereas brown (grey) lines show the symmetrically evaluated 
g^^\—z/2, z/2) ofEq. (9). Plotted graphs correspond to tem- 
peratures r/r^ « (i) 0.1, (h) 0.5, (hi) 0.9, and (iv) 2.8, with z 
scaled throughout to the temperature-dependent system size 
Rtf{T). (b) Value of the normalised first-order correlation 
function evaluated at 2; = Rtp (T) /2 from the trap centre, as a 
function of scaled temperature T/T^, obtained from (j'^'(0, 2) 
(brown/grey circles), and g'^' (—2/2, z/2) (black squares), (c) 
Dependence of Rtf{T) on scaled temperature T/T^. Plotted 
lines in (b)-(c) are a guide to the eye. 
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FIG. 3: (color online) (a) (i)-(iv) Symmetrically computed 
spatial correlation function g'^' (—2/2, 2/2), where 2 is scaled 
to the temperature-dependent Thomas- Fermi radius Rtf(T) 
(solid brown/grey), versus numerical fit by the function /(z) 
of Eq. (10) (dashed black), for the graphs of Fig. 2(a). (b) 
Temperature dependence of the coherence length I/J-oh (black 
circles) obtained by above fit, with L'^ati scaled to Rtf{T), and 
temperature to T^. The dashed brown (grey) line highlights 
a coherence equal to the system size, i.e. L'^^^^ = i?Tp(r). 
(c) Corresponding dependence of the 'crossover' parameter ( 
(black circles) determining the relative importance of expo- 
nential and gaussian contributions to the spatial correlation. 
The dashed brown (grey) line separates regions of predomi- 
nantly gaussian behaviour ((" 2> 1), to those of exponential 
behaviour. Plotted solid lines in (b)-(c) are a guide to the 
eye. 



C. Exponential vs. Gaussian Profiles 

The presented analysis thus far has focused only on 
the 'global' loss of coherence with increasing tempera- 
ture, without paying attention to the precise shape of 
the correlation function. The latter is known to vary 
from exponentional to gaussian with decreasing tempera- 
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ture, indicating the crossover from quasi-condensation to 
'true' condensation |2^. In this section we study the tem- 
perature crossover in the shape of this function in more 
detail, by means of our previously computed symmetric 
correlation function g^^' (— z/2, z/2). Since g*-^^ is known 
to be gaussian at T = and exponential at temperatures 
T > T^, we fit it, in the relevant region z < 2R'tf{T), 
by a function f(z) which provides a smooth crossover 
between the two shapes, of the form 



(10) 



The parameter ^ appearing above is a measure of how 
'gaussian' or 'exponential' a particular profile is, while 
the chosen fit provides a unique definition of the coher- 
ence length, L'^^^, which decreases monotonically with 
increasing temperature, irrespective of the change in the 
profile's shape. 

The correlation functions z/2, z/2) for the four 

different temperatures discussed earlier are shown by the 
brown (grey) lines in Fig. 3(a); the corresponding fits 
based on Eq. (10) are plotted on the same figure by 
dashed black lines. Such fits enable us to study the de- 
pendence of the 'crossover parameter', C, and the coher- 
ence length, -/j^oh' temperature. In particular, the 
scaled coherence length L'^^^^/ Rtf{T) and the 'crossover' 
parameter C are plotted against scaled temperature T/Trp 
in Fig. 3(b)-(c). Note that a value of C = corresponds 
to a purely exponential correlation function, with the 
gaussian limit reached for C ^ 1. 

At low temperatures, T < T^, we find that the co- 
herence length ipoh comparable to, or larger than the 
quasi-condensate spatial extent, R^p{T), for that par- 
ticular temperature, signalling the appearance of 'true' 
condensation. In this regime, the gaussian contribution 
has a larger weight, than the exponential term. More- 
over, we find that the ratio in their relative contributions 
decreases with increasing temperature, with ^ « 1 at the 
'crossover region'. Finally, as the temperature increases 
further, the system enters the quasi-condensate regime 
and the correlation function essentially acquires an expo- 
nential profile for temperatures T > T^, consistent with 
experimental observations. 



IV. HIGHER ORDER (DENSITY-DENSITY) 
CORRELATION FUNCTIONS 



In this section we consider higher order correlation 
functions, which contain information about density- 
density correlations. These are considered at fixed 
points z and at equilibrium, with the subsequent anal- 
ysis focusing on the correlation functions ^'^^^(z) = 
g^'^'>{z, z, z, z;teq), and g^^\z) = g'^^'' (z, z, z, z, z, z;teq) of 
Eqs. (6)-(7) as a function of z for various temperatures. 



A. Dependence on Temperature 

The behaviour of g''^\z) and g^^\z) against z is shown 
for various temperatures in Figs. 4(a) and (b), with po- 
sition scaled to the zero-temperature Thomas-Fermi ra- 
dius, i?TF(0). As expected, these functions have a lower 
bound of 1, occuring deep within the quasi-condensate 
region (z ^ i?TF(0)), and at sufficiently low tempera- 
tures, demonstrating full coherence. The corresponding 
upper bounds for these functions are ^'■^-'(O) = 2! = 2, 
and g*-'^^(0) = 3! = 6, consistent with the observed values 
in 3D BECs US |M El El- These values correspond 
to complete absence of coherence, and such incoherent 
thermal atoms are typically located outside the quasi- 
condensate region. 

In particular, as T — > 0, both functions g^^^'iz) and 
g^^^z) tend towards step-like functions g^"^z) = 1 + 
(n! — l)6(z) centered around i?TF(0), with Q{z) — for 
z < i?TF(0) and 1 for z > i?TF(0). As the temperature 
increases, we observe the following features, which are 
consistent with the results of Firstly, the central 

values g*-^^(0) and g'^'^\0) increase due to the increas- 
ing thermal component located at the trap centre; their 
dependence on temperature is shown in Fig. 4(c). Fur- 
thermore, the crossover between the quasi-condensate- 
dominated region z < Rj^piT), and the purely thermal 
region z > Rtf{T) becomes smoother, due to the in- 
creased presence of thermal atoms over the entire trap 
extent. Finally, the location of this crossover is shifted 
to smaller values of z/i?TF(0), consistent with the de- 
crease in the spatial extent Rtf{T) < i?TF(0) of the 
quasi-condensate with increasing temperature. 

Our findings for the temperature dependence of g^^\0) 
in Fig. 4(c) are further compared to the analytical result 
of Kheruntsyan et al. pa . l57l| . In their work, the interac- 
tion strength is parametrised in terms of the parameter 
7 = mg/h^n, where m is the atomic mass, g the one- 
dimensional coupling constant, and n the density of the 
gas. In the regime 7 <C 1 considered here, they find the 
following temperature dependence (see Eq. (5.10) in Ref. 



3(^^(0) -1 + 



4%/2 



(11) 



where is the degeneracy temperature, defined by 
Td — Nhwz, N is the number of quasi-condensate atoms 
and oJz is the confining frequency. We apply this an- 
alytical formula to our numerical results by extracting 
the relevant parameter (T /Td) from our stochastic sim- 
ulations. The corresponding temperature dependence of 
g^^-* (0) predicted by Eq. (11) is thus plotted by the brown 
(grey) triangles in Fig. 4(c). Thus there is good agree- 
ment between the stochastic and analytical theories in 
the temperature range considered here. However, the 
stochastic theory predicts an additional slight reduction 
in coherence, as manifested by the higher value of g^"^^ (0), 
with the extent of this reduction increasing with increas- 
ing temperature. 
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FIG. 4: (color online) (a)-(b) Spatial dependence of the 
(a) second, g^^\z), and (b) third order, g^'^\z), equilibrium 
spatial correlation functions at the same point z, for var- 
ious temperatures; from bottom to top T/T^ ~: (i) 0.1, 
(ii) 0.7, (iii) 1.7, and (iv) 3.5. Distances are scaled to 
the zero temperature system size _Rtp(0). (c) Correspond- 
ing values of ^^^•'(O) = g*^' (O, 0, 0, O; te,) (black circles) and 
^'^'(O) = p'^'(0, 0,0,0, 0,0; feq) (black squares) at trap centre 
as a function of temperature. Brown (grey) triangles indicate 
corresponding results for g^^\0) based on Eq. (11). Plotted 
lines in (c) are a guide to the eye. 



B. Dependence on Interaction Strength 

Next, we investigate the dependence of the three-lowest 
order correlation functions on the effective ID interac- 
tion strength at fixed temperature, with our analysis 
again restricted to the weakly-interacting regime 7^1. 
As mentioned earlier, the ID interaction strength can 
be parametrised in terms of the ID coupling constant 
K = (m/27r?i)a3DW_L. This parameter can be changed 
either by tuningthc 3D scattering length via Feshbach 
resonances '66', , or by modifying the transverse con- 
finement, while keeping all other parameters fixed. Since 
an increase in the interaction strength causes a reduction 
of coherence, as demonstrated explicitly below, we have 
chosen to discuss here a rather low temperature example, 
with T = O.ST^, for which there is appreciable coherence 
in the system. 

The dependence of the three lowest order correlation 
functions on coupling constant k is shown in Fig. 5, 
with the coupling constant used thus far in our discus- 
sion, and in Figs. 1-4, henceforth denoted by kq. Since 
temperature is kept fixed for all curves shown in this 
section, we have chosen here to scale distances to the 
quasi-condensate spatial extent at the given temperature 
T = O.ST^, based on the coupling constant kq. For this 
particular coupling constant, the corresponding first or- 
der correlation functions were shown in Figs. 2(a)(ii) and 
3(a)(ii) (note the different plot ranges), whereas higher 



FIG. 5: (color online) Dependence of the three lowest order 
correlation functions on the effective ID coupling constant k 
at equilibrium for a fixed low temperature T ~ 0.5T^: (a) 
Dependence of g^^\—z/2, z/2) on «: as a function of position, 
for (from top to bottom) k/kq — 1, 2, 3, 5 and 10, where kq is 
the value of the effective ID coupling constant used in Figs. 
1-4. (b) Dependence of g^^\-z/2, z/2) at z = 7?TF(r)/2 
on k/ko, with the dashed line simply connecting the data 
points, (c)-(d) Second and third order spatial correlation 
functions at the same point and at equilibrium as a function 
of position, for an interaction strength (from bottom to top) 
k/ko = 1,5, 10. Here, the value of the ID coupling constant 
K has been changed by modifying the transverse confinement 
from ojx = 27r X 120Hz, to the values u}'^/2n = 240, 360, 600, 
and 1200 Hz. In this figure, all positions have been scaled to 
the spatial extent of the system at T = 0.5r,^ for a coupling 
constant k = ko- 



order correlation functions were plotted by the dashed 
black lines in Figs. 4(a)-(b). 

The change in the correlation function g''-^'' {—z/2, z/2) 
as a function of position for different values of the ID 
coupling constant is shown in Fig. 5(a). As evident, an 
increase in k, from its previously considered value kqj 
leads to a decrease in the coherence of the system. This 
is clearly portrayed in Fig. 5(b) showing the value of the 
correlation function at z = Rt-p(T)/2 as a function of 
K. This curve further demonstrates that full coherence 
is asymptotically approached in the limit of negligible 
interactions, i.e. as k — > 0, as anticipated. 

Figs. 5(c)-(d) show corresponding dependence of 
5^^-' (z) and g'-^^ (z) on z for different interaction strengths, 
with top curves corresponding to a ten-fold increase in k. 
Our analysis reveals that an increase of the interaction 
strength leads to an increase in the central values of the 
correlation function g^^^^ (0) and g^^^^ (0), and a more grad- 
ual crossover to incoherent behaviour outside the quasi- 
condensate region. Thus, qualitatively, an increase in the 
interaction strength produces the same effect on the sys- 
tem coherence as an increase in the temperature of the 
system. 
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V. DIRECT DETERMINATION OF 
QUASI-CONDENSATE PROFILES 

The calculation of density-density correlations dis- 
cussed previously enables the precise determination of 
the quasi-condensate density, to which we now turn our 
attention. 



A. Methodology 

To be able to identify the quasi-condensate part in our 
density profiles, we use the fact that the Langevin field 
<i>(z,t) essentially contains the physics of the Bose field 
operator ^{z,t). Such a correspondence has already been 
used in numerically evaluating the three lowest order cor- 
relation functions g^^\- • • ), f/^^H' ' ' )i d'^^H' ' ' ) from 
Eqs. (5)-(7) in Sees. Ill and IV. At a given time, the 
total atomic density is obtained via 

n(z) ^ ^ , (12) 

where the latter averaging is performed over different re- 
alizations of the noise. 

In order to separate out quasi-condensate and thermal 
contributions, we impose the usual 'decomposition' of the 
Bose field operator into a 'mean field part' ipo{z) and a 
'fluctuating part' S{z), i.e., 

= + 5{z) = M^^) + k^) ■ (13) 

This directly enables us to split the atomic profile of Eq. 
(12) into a 'quasi-condensate part', nQc(^), and a 'ther- 
mal part', nT{z), via 

n(z) = nQc(z) + rniz) = \M^)f + (<5t(z)5(z)) . (14) 

We now seek to re- write these two components in terms 
of expressions involving same-time averages of multiple 
Bose field operators, since the latter quantities can be 
routinely evaluated within our numerical scheme. In par- 
ticular, we note that the average over four Bose field op- 
erators at the same spatial coordinate (and same time) 
can be written as 

{i'\z)i'\z)i'{z)i'{z)) = \Mz)f 

+4 \M^)f {6\z)6{z)) + 2{5Hz)S{z)f (15) 

where we have used Wick's theorem in the form 

{S\z)S\z)Siz)S{z)) = 2(S\z)S{z))^ . (16) 

At the same time, we make use of the following corre- 
spondence 

(^'t(z)*t(z)*(z)#(z)) ^ ($*(z)$*(z)$(z)$(z)) 

= iimi") (17) 



So, from Eqs. (12)-(17) it is easy to see that, in the region 
z < i?TF(T) where a quasi-condensate exists, and there- 
fore the above splitting into mean field and fluctuations 
is meaningful, the following relation holds 

\Mz)\'^ = 2(*^(z)*(z))2 - {ii\z)i'\z)i'{z)i'{z)) . 

(18) 

The quasi-condensate density nQc(^) = |V'o(-z)|^ can 
hence be accurately determined in the range z < i?TF (T) 
from our numerical scheme via 

UQciz) = y^2(*t(z)*(z))2 - {^^z)i'Hz)^{z)i'iz)) 
- v/2(|<l>(z)P)2-(|$(^)|4). (19) 

Correspondingly, the profile of the thermal cloud, 
nT{z) = {S''{z)S{z)) is directly obtained via 

nriz) ^ {mz)f) - UQciz) . (20) 

This approach has already been used (but not com- 
mented upon) in parallel work discussing quasi- 
condensate growth on an atom chip '58|. 

Note that an alternative approach to determine quasi- 
condensate density profiles has been recently discussed 
in l45l . based on the projected Gross-Pitaevskii scheme 



(a) (b) 




(c) (d) 




1 2 3 1 2 3 



Z / Rtf(T) 

FIG. 6; (a)-(b) Total density profiles, scaled to the T = 
peak central density, riff'', as obtained from the stochastic 
approach (black), and mean field theory with density fluctu- 
ations (brown/grey) corresponding to the correlation func- 
tions plotted Fig. l(a)-(b). (c)-(d) Corresponding quasi- 
condensate (dashed) versus thermal cloud profiles (solid). 
Insets: Thermal cloud profiles around their maxima, z — 
Rtf (T) . All graphs are scaled to the temperature-dependent 
Thomas- Fermi radius _Rtf (T) for the particular temperature, 
and are based on the coupling constant kq of Figs. 1-4. 
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B. Results 

To apply the above method, we start by plotting the 
total density profiles at two different temperatures, above 
and below in Fig. 6(a)-(b). These figures depict the 
good agreement between the stochastic and mean field 
approaches commented upon in earlier work |23j | . In this 
paper such an analogy is extended much further, by di- 
rectly comparing quasi-condensate and thermal profiles 
predicted by these two theories, as shown in Fig. 6(c)-(d). 
The quasi-condensate density is evaluated in the stochas- 
tic theory by Eq. (19) within the range < ^ < Rtf{T), 
and is zero elsewhere. 

In order to make a meaningful comparison between 
these distinct approaches, at any given temperature, 
the same temperature-dependent Thomas-Fermi radius 
is used for both theories, with Rtf{T) evaluated from the 
mean field theory. Dashed lines in Figs. 6(c)-(d) depict 
the quasi-condensate, as evaluated from the stochastic 
(black) and the mean field theory (brown/grey), whereas 
solid lines show corresponding thermal cloud distribu- 
tions. The latter are highlighted in the insets around 
their respective maxima, located at the quasi-condensate 
edges z = Rtf{T). The agreement between the pre- 
dicted profiles is excellent, and remains so over the entire 
temperature range considered. 

The approach presented here further enables the direct 
determination of the equilibrium quasi-condensate frac- 
tion at any temperature. In particular, we note that the 
quasi-condensate fraction is found to decrease approxi- 
mately linearly with increasing temperature. 

C. Application to Experiments 

Experiments with quasi-condensates typically rely on 
bimodal density fits to determine important parameters 
such as the quasi-condensate spatial extent and quasi- 
condensate fraction at a given temperature. The proce- 
dure described above provides a direct alternative way of 
obtaining an unequivocal, theory-independent determi- 
nation of quasi-condensate and thermal density profiles 
in quasi-lD experiments. 

In order to perform such an analysis experimentally, 
one must simultaneously obtain both the total atomic 
density profiles n{z) — (|<I>(z)p), as done routinely, 
as well as the density-density correlations {n^{z)) — 
(|$(2:)|^) at various points across the entire trap. This 
experiment should ideally be performed in situ, which is 
within current experimental reach [s^, in order to avoid 
the coupling of density and phase fluctuations occuring 
in the usual time-of- flight expansion stage 0, |^ . 

VI. CONCLUSIONS 

In conclusion, we carried out a systematic analysis 
of the dependence of the three lowest order correlation 



functions of a one-dimensional ultracold atomic gas on 
temperature and effective interaction strength, by means 
of a stochastic Langevin approach in the 'classical' ap- 
proximation. Our discussion was limited to the weakly- 
interacting regime (7 ^ 1), which features a smooth 
temperature crossover from an incoherent gas, to a quasi- 
condensate, and, finally, at sufficiently low temperatures, 
to a 'true' condensate. This work complements earlier 
work in this area, which studied the interplay between 
density and phase fluctuations. 

Results for the off-diagonal normalized first order cor- 
relation function were shown to be practically indistin- 
guishable from those of a mean field theory which is valid 
in low dimensions and explicitly includes density fluc- 
tuations. This should be contrasted to theories which 
a priori ignore density fluctuations, with such theories 
shown to predict a larger amount of coherence than is 
physically present. Two different approximations for the 
evaluation of this correlation function were used, corre- 
sponding to different interference experiments that could 
be performed, and their profiles were shown to be very 
similar within the appropriate spatial range. We further 
investigated the temperature crossover in the shape of 
the correlation function, which evolves from gaussian to 
exponential with increasing temperature, and used this 
to obtain a monotonically varying coherence length. 

Moreover, second and third order correlation functions 
at the same point were investigated as a function of tem- 
perature and interaction strength. The observed decrease 
in coherence when either of these parameters is increased 
was shown to be consistent with other results in the liter- 
ature. Importantly, consideration of density-density cor- 
relations was shown to lead to an accurate determina- 
tion of the quasi-condensate and thermal density profiles. 
This approach could be used in experiments, to provide 
a more accurate determination of quasi-condensate pro- 
files and fractions, offering an alternative to the bimodal 
fitting schemes currently used. 

In addition to providing valuable information for equi- 
librium profiles, to which the present analysis was re- 
stricted, the stochastic approach is well-suited for dis- 
cussing growth of coherence in quasi-one-dimensional 
Bose gases, an issue that will be addressed in subsequent 
work. 
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APPENDIX A: LOW-DIMENSIONAL MEAN 
FIELD THEORIES 

This Appendix summarizes the modified finite tem- 
perature mean field theory of Andersen, Al Khawaja and 
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Stoof [231 ■ This theory extends alternatiye approaches 
suited to low-dimensional gases |^ by proyiding an 
ab initio self-consistent treatment of density and phase 
fluctuations, which yields an equation of state free of both 
infrared and ultrayiolet diyergences. This theory is yalid 
both for homogeneous and trapped systems in all dimen- 
sions, and has been shown to be consistent with well- 
known results in the appropriate limits 23]. Further- 
more, such an approach enables a direct determination 
of quasi-condensate density profiles. 

In the harmonic trap, one obtains an appropriate 
finite temperature nonlinear Schrodinger equation for 
the quasi-condensate, coupled to the Bogoliubov-de- 
Gennes equations for the excitations. Solying these 
self-consistently within the Thomas-Fermi approxima- 
tion yields a mean field approximation for the quasi- 
condensate spatial extent. This is referred to as the 
temperature-dependent Thomas-Fermi radius, and is de- 
noted by RtfIT). This quantity depends on the quasi- 
condensate depletion n'{z) — nT{z), via R^p{T) = 
■y/2/i'/r7ia;^ where ^' = /i — 2Kn'(0) is the 'renormalizcd' 
chemical potential. Here k is the effective ID coupling 
constant which is fixed by the atom mass, m, the 3D 
scattering length, osd, and the transverse co nfining fr e- 
quency u!± via k = a3D/2TT£']_, where l± — ^/h/muTi is 
the harmonic oscillator length in the transverse direction 
|65|. Density profiles are obtained from the above equa- 
tions in the local density approximation, with thermal 
profiles outside the quasi-condensate region obtained by 
solving the equation of state for the normal gas. 

The first-order correlation function can be written as 
5(1) (0, z) = exp (-([x(^) - x(0)]')/2) , where x{z) is the 
operator for the phase, defined by expressing the Bose 
field operator vE'(z) in the density-phase representation 
as ^(z) = y/ n(z)expU'y(z)}. In the theory of Ander- 
sen et al. [23, where density fluctuations are explicitely 
taken into consideration, the above exponent assumes. 



in the purely one-dimensional geometry considered here, 
the following form 



(Al) 



/ P,{z/RMT)) ^ 



Pj[z) are 
J, with A 



iViT 



Legendre polynomials of or- 
. y/{j + l/2)fi'/hiuj, and 

l/2)hujj/ fi')/2. The frequencies are 
is the harmonic 



Here 
der 
B, - 

given by ujj = + l)/2 oJz, where 

oscillator length corresponding to a longitudinal confin- 
ing frequency lo^- N{hujj) is the usual Bose distribution 
function. 

The appearance of the temperature-dependent quasi- 
condensate size i?TF (T) in the prefactor ensures that den- 
sity fiuctuations are explicitly maintained in the above 
expression. The general expression quoted above can be 
readily reduced to the conventional theory which ignores 
quasi-condensate depletion 0, by replacing _Rtf (T) 
by the corresponding zero temperature quasi-condensate 
size R-jp{0) at the same total atom number, and ignoring 
the Bj contributions. Since R'tp{T) cx y/Jj/, the former 
step is equivalent to replacing the renormalizcd chem- 
ical potential at temperature T by the corresponding 
zero-temperature one for the same total atom number. 
Then, the 'classical' approximation N(hu!j) « kBT/hwj 
leads to the definition of the characteristic temperature 

For more details on the implementation of this theory 
to trapped gases, the reader is referred to [2^ I23ll2^ I2,'t| . 
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